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Abstract In some previous articles, we defined several partitions of the total kinetic energy 
T of a system of classical particles in into components corresponding to various modes 
of motion. In the present paper, we propose formulas for the mean values of these compo- 
nents in the normalization T = \ (for any d and A') under the assumption that the masses of 
all the particles are equal. These formulas are proven at the "physical level" of rigor and nu- 
merically confirmed for planar systems (d = 2) at 3 $J A' 100. The case where the masses 
of the particles are chosen at random is also considered. The paper complements our article 
of 2008 [Russian J Phys Chem B 2(6):947-963] where similar numerical experiments were 
carried out for spatial systems (d = 3) sA3 '^N ^ 100. 
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1 Introduction 

The integral characteristic of motion in a system of classical particles is the total kinetic 
energy T . However, for quite different types of motion, the value of T can obviously be 
the same. For instance, consider a system of two particles with the fixed center-of-mass 
and let r = r{t) be the vector connecting the particles. Depending on the forces and initial 
conditions, this two-particle system can vibrate (if the direction of r is not changing) or 
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rotate (if the length of r is not changing), and the kinetic energy in both the cases may 
attain any positive value. However, these two kinds of motion are in fact extremes: a typical 
motion is a mixture of vibrations and rotations. Of course, in the case of two particles, the 
total kinetic energy T can be straightforwardly and naturally decomposed as the sum of two 
terms corresponding to vibrations and rotations. Indeed, let r = dr/dt = r|| +rx where 
is the component parallel to r and Fx is the component orthogonal to r. Then 

rr _ m\m2 ij.|2 _ jVib _|_ y,i-ot 



2{rn\ +WI2) 

where 



= ^ ' ^ J i-ii = rcos^e 

2(mi +1112) 



is the vibrational energy and 



rl"Ot 



m\m2 



2(mi +m2] 



r 



= Tsin^e 



is the rotational energy (here mi, mi are the masses of the particles and Q is the angle 
between r and r). 

In the case of three or more particles, there are much more kinds of motion: one has 
to distinguish rotations of the system in question (the ''cluster" or ''aggregate") as a whole, 
changes in the principal moments of inertia, various rearrangements of particles in the clus- 
ter, etc. In this situation, it becomes a rather non-trivial and ambiguous task to define the 
components of T corresponding to various modes of the motion. This problem has been 
discussed in the literature for decades and several approaches to kinetic energy partitioning 
have been proposed; see e.g. the well known papers [1,2 3, 4|, the recent studies by Marsden 
and coworkers El|6], and references therein. 

The importance of exploring the contributions of various motion modes to the total ki- 
netic energy T stems, to a large extent, from the fact that such energy components may most 
probably be used as effective global indicators of dynamical features and critical phenomena 
(e.g. phase transitions) in classical clusters. If these energy terms are defined in a sufficiently 
"symmetric" and "invariant" manner and can be computed automatically and fast from the 
coordinates and velocities of the particles, then one expects to be able to straightforwardly 
detect structural metamorphoses in the cluster by observing abrupt changes in the way the 
total kinetic energy is distributed among the modes. In the case of large clusters, this is cru- 
cial for applications because it is much easier to trace a few indicators than to analyze a huge 
collection of data pertaining to all the particles. 

Starting in 2002, we have published a series of papers ll7]|8ll9] |T0l[m where we pro- 
posed and preliminarily tested a number of such global indicators on the basis of the so- 
called hyperspherical approach to cluster dynamics. In particular, several partitions of the 
total kinetic energy T with various amazing features were defined in [9 | for particles in the 
conventional three-dimensional space and in [ lOJ for particles in the Euclidean space of an 
arbitrary dimension d. Moreover, our long paper IIIOI contains also rigorous mathematical 
proofs of many properties (including the invariance under certain group actions) of the terms 
of these partitions. A refinement of one of the partitions was performed in 1,11] and four new 
terms were introduced. Some perspectives of the partitions in question are discussed from 
the general viewpoint of the methods of molecular dynamics (with the particles being atoms 
or ions) in the short reviews [12, 13 14|. Applications include studies of such phenomena 
and processes in chemical physics as phase transitions in small neutral argon clusters Arn 
with n = 3, 13, 38, 55 ll8l ll5lfT6i , dynamics and thermodynamics of small ionic argon clusters 
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Ar+ with n = 3, 6, 9 1 17], the prototypical exchange reaction F + H2 — > H + HF of atomic 
fluorine and molecular hydrogen [15], and ultrafast relaxation dynamics of krypton atomic 
matrices doped with a nitrogen monoxide molecule [18 19 1 (some of these applications are 
surveyed in |14|). These studies have confirmed the usefulness of the kinetic energy parti- 
tions for examining dynamics of classical nanoaggregates. 

However, neither the theory developed in the articles f8^'97'TOl nor the applications con- 
sidered in f8lll5lll6lfr7lll8lll9l allow one to conclude how the energy partitions defined in 
I MTOl look like for "typical" systems (roughly speaking, for a random choice of the coordi- 
nates and velocities of the particles) or how the statistics of the terms depends on the number 
A' of the particles and on their masses. The statistics of the kinetic energy components for 
systems in the conventional space M.^ was studied numerically in our paper 0111 in the range 
3 A' ^ 100 for two extreme situations: in the case where all the particles have equal masses 
and in the case where the masses are chosen at random. The main observation of Hill is that 
in the situation of equal masses, the mean values of almost all the terms of the partitions in 
the normalization J = 1 are very simple (in fact, linear fractional for N ^ 4) functions of A'^. 
For some of the terms, the paper iHT l proposed also generalizations of the formulas for the 
mean values to an arbitrary dimension d of the space. 

The aim of the present article is fourfold. First, we rigorously prove the properties of the 
new energy terms S™'' ^1"' ^2" announced in our previous paper [iTl. Second, in the 
situation of equal masses, we suggest formulas for the mean values of all the energy com- 
ponents (except for the so-called unbounded ones) for A' particles in with an arbitrary d: 
the mean values turn out to be simple rational functions of d, N, and mm{d,N — 1). Third, 
we prove these formulas at the "physical level" of rigor. Fourth, we examine numerically the 
statistics of all the energy terms, both in the situation of equal masses and in the situation of 
random masses, for A' particles on the plane in the same range 3 ^ A' $S 100 as in 1111 . 
In the situation of equal masses, our numerical experiments do confirm the formulas for the 
mean values of the energy terms. 

It is worthwhile to note that two-dimensional physics is a well developed and rapidly 
progressing field of science, see e.g. the works II20II21II22]|23II24||251|26II27I I28''291 and ref- 
erences therein (by the way, the preprint [24] contains 680 references). Of course, here the 
words "two-dimensional" can be used either in the mathematical sense or in the sense of 
monolayer films. In particular, the Nobel Prize in Physics of 2010 was awarded jointly to 
Andre K. Geim and Konstantin S. Novoselov "for groundbreaking experiments regarding 
the two-dimensional material graphene" Il26ll27il30i |. Moreover, a paradigm was recently 
proposed in which the actual space- time is a fundamentally (1 -|- 1) -dimensional universe 
but it is "wrapped up" in such a way that it appears higher dimensional (say, 2-1-1 and 3-1-1) 
at larger distances II31II32II33I . 

The paper is organized as follows. Some fundamental theoretical concepts related to our 
approach to kinetic energy partitions are recalled in Sections [2] and [3] while the partitions 
themselves are defined in Section |4] Almost all the material of Sections |2}^ is in fact con- 
tained in our articles [9 10] and is reproduced here for the reader's convenience. The new 
energy terms E™, E2 of 1 1 1 1 are defined and explored in Section|5] The formulas 

for the mean values of the bounded energy components for A^ particles of equal masses in 

are proposed in Section|6]and substantiated in Section|7] The numerical experiments for 
d = 2 and their results are described in Section|8] Concluding remarks follow in Section|9] 
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2 Instantaneous phase-space invariants 

Consider a system of A' classical particles in the Euclidean space K.'' {d,N €N= {l;2; ...}) 
with masses m\,m2, . . .,mi^. Let ri = ri (f), r2 = r2(f), . . . , Fat = Ti^{t) be the radii vectors 
of these particles with respect to the origin. Introduce the notation 

N 

M=Y^ Ma, qa = {ma/Myl\a 

a=\ 

(\ ^ a ^ N). In terms of the total mass M of the system and the mass-scaled radii vectors 
qa, the total kinetic energy T and the total angular momentum / of this system are expressed 
in an especially simple way: 

\ N N 

r=^Ema|r„p = ^ j:iq„P, (1) 

N N 
= L ■^O' -^'V = L '«a(r,;ar;;a -r^;„r,;a) =MY^ {<li;a<\r,a-<\r,a<\i;a) (2) 

l<ii<j^d a=l a=l 

(1 /' < J ^ d). Here r,;c( and qi-a denote respectively the ;th components of the vectors 
Ta and qa and, as usual, the dot over a letter means the time derivative. For d = 3, the 
formulas ([T} and (|2]i give the conventional total kinetic energy and total angular momentum 
of a system of A' particles. Thus, the mass-scaled radii vectors q^ (f) provide a more adequate 
description of the current state of a system of classical particles than the radii vectors (t) 
themselves. 

Definition 1 The position matrix of a system of A' classical particles in R'' is the d x N 
matrix Z whose columns are the mass-scaled radii vectors qi , q2, . . . , qN- 

Of course, the position matrix at any given time instant t depends on the Cartesian co- 
ordinate frame chosen. The choice of another coordinate frame (with the same origin) is 
described by a transformation Z(f) ^ RZ{t) with an orthogonal d x d matrix R 6 0{d). 
On the other hand, according to the general duality concept for the physical space and the 
abstract "kinematic space" iQl llOllllI (which are respectively M'' and M.^ in our case), one 
can also consider transformations of the form Z(f) Z{t)Q* with orthogonal N x N matri- 
ces Q 6 O(A') where the asterisk designates transposing. Such transformations correspond 
to changes in the "type" of the coordinate frame. For instance, the passage from Cartesian 
coordinates to Jacobi or Radau-Smith coordinate frames (see e.g. the papers P1 ll0ll34ll35l 
and references therein) is equivalent to a multiplication of the position matrix from the right 
by a suitable fixed orthogonal matrix [ 10 1. 

Example 1 Suppose that the A' ^ 2 particles in question move without external forces and 
their center-of-mass coincides with the origin. Then in any Cartesian coordinate frame 

N N 

m]l'^qa=M'^l^ ^ mar a =0. 
a=l a=l 

Choose an arbitrary matrix Q 6 O(A') with the entries in the last row equal to 

QNa = {malM)^l^, 1 ^ a ^ ^. 

Then the last column of the matrix Z{t)Q* is zero at any time moment t (and one may 
therefore regard the "kinematic space" to be M'*'^'). 
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We are interested in characterizing systems of classical particles in by various quan- 
tities that are determined, at any time instant to, by the positions Z(fo) and velocities Z{to) 
of the particle^ at this time moment to onl}0 (rather than by the whole trajectory {Z(f) | 
'begin ^ ' 'end}) and are invariant under orthogonal coordinate transformation^ both in the 
physical space and in the "kinematic space". Being inspired by the discussion above and, in 
particular, by Example [T] one arrives at the following mathematical model. 

Let rf,n £ N. Consider the space (R"'^")^ of the pairs of real dxn matrices. We will 
write down these pairs as {Z,Z). On the space (W^^")^, there acts the group 0{d) x 0(«): 

{R,Q){Z,Z) = {RZQ*,RZQ*), ReO{d), QeO{n). (3) 

Definition 2 An instantaneous phase-space invariant of systems of classical particles in R"' 
is a collection of functions 

/rf,„ : Mrf,„ X R''^'" ^ R, Md.„ C n > no 

possessing the following two properties. 

First, the definition domains Mrf „ x R'''^" of the functions /j „ and the functions „ 
themselves are invariant under the action (|3} of the group 0{d) x 0(n): 

zeMa.n =^ RZQ*eMd,, & fd.„{RZQ*,RZQ*)=fd,,{z,z) 

for any Z 6 R''''" (« ^ «o), R 6 0(J), and Q 6 0(w). 

Second, the value of the function remains the same if one augments both the matri- 
ces Z and Z by the (« + l)th column equal to zero: 

zeMrf,„ =^ (Z0)eMrf,„+i & /^,„+i((zo),(zo)) =/rf,„(z,z) 

for any Z e R"''"' (« > no). 

Most of the instantaneous phase-space invariants we will deal with will depend on the 
parameter M (the total mass of the system of particles). 

Example 2 Let the center-of-mass of a system of N ^ 2 particles in R'' coincide with the 
origin. As was explained in Example[T] one can choose a (non-Cartesian) coordinate frame ^ 
in which the last column of the position matrix Z(f ) of such a system is identically zero. The 
first A' — 1 columns of Z(t) constitute the d x {N — I) matrix z''^''™'''' (t). Any instantaneous 
phase-space invariant / of the system in question can be calculated either using the d x N 
position matrices in the initial Cartesian coordinate frame (n = A') or using Ihe d x [N — I) 
reduced position matrices z"^duced jj^g coordinate frame ^ {n = N — I). 

In the sequel, it will be convenient to introduce the notation 

m = mm(d,n). 

As is very well known (see e.g. the manuals Il36ll37ll38ll39l ), any dxn matrix Z 6 R''^" can 
be decomposed as the product of three matrices 

Z = DTX*, DeO{d), XeO{n), (4) 

' Whence the words "phase-space" in Definition[2]below. 
^ Whence the word "instantaneous" in Definition|2]beIow. 
^ Whence the word "invariant" in Definition [2] below. 
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where all the entries of the J x « matrix T are zeroes with the possible exception of the 
diagonal entries: 

^11 = ^l, ^22 = i§2, • • • , ^Tim = ^m, <§! !^ "§2 ^ ' ' ' ^ "^m ^ 0. 

The representation © is called the singular value decomposition (SVD) of the matrix Z. 
The numbers , (§2, • ■ • , are called the singular values of the matrix Z and are determined 
uniquely by Z although the orthogonal factors D and X in the equality ([4]i are not. If n ^ d 
then the d singular values of Z are the square roots of the eigenvalues of the symmetric dxd 
matrix ZZ*. If n ^ d then the n singular values of Z are the square roots of the eigenvalues 
of the symmetric nx n matrix Z*Z. 

It is clear that each singular value of Z is an instantaneous phase-space invariant; to be 
more precise, for each I ^ i ^ d, the ith singular value : Kf^^" — >• R is an instantaneous 
phase-space invariant for n ^ /. More generally, any function of the singular values of Z 
is an instantaneous phase-space invariant, and conversely, any instantaneous phase-space 
invariant independent of Z is a function of the singular values of Z. 

Along with the SVD of a matrix Z, it is often expedient to consider a signed singu- 
lar value decomposition (signed SVD) |10,40|. A signed SVD of a d x n matrix Z is a 
representation (|4) where again all the entries of the d x n matrix Y are zeroes with the pos- 
sible exception of the diagonal entries Tii,Y22, ■ ■ ■ ,Tmm^ but it is no longer assumed that 
Yii ^ Y22 ^ ■ ■ ■ ^ Txnm ^ 0. Of course, the numbers iTn |, |r22|, • • • , l^lnml constitute in this 
case an unordered collection of the singular values of Z. The differentiability properties of a 
signed SVD are usually better than those of the SVD (see the papers II10II40I and references 
therein). 

Recall that the standard Frobenius inner product on the space R''^" of real dxn matrices 
is defined by the formula 

{Z\Z'') = Tr[Z«(Z^)*] = t^f^Zf^' (5) 

a=l 1=1 

where Tr denotes the trace of a square matrix. The corresponding matrix norm || || given by 

n d m 

||zf = (z,z) = Tr(zr) = l^zl =Y^e, 

a=\i=\ (7=1 

is called the Frobenius norm Il36ll37l[39l , the Euclidean norm II37II38I . the l2-norm 1371 . the 
Schur norm |37|, the Hilbert-Schmidt norm (37], or the spherical norm f381. 

In the context of the pairs (Z,Z) £ (R''^")^, the instantaneous phase-space invariant 

P = I|Z|| 

is called the hyperradius (of the system of particles whose position matrix is Z). 

Remark 1 In the sequel, we will use the following notation. Given a pair (Z,Z) £ (R''^")^ 
and a certain object (number, vector, etc.) X = X{Z) dependent on Z, we will define X. = 

X{Z,Z) as 

*=^^(3(0) , 

f=0 

where 3 : {-e,e) ^R''"" (0 < e «: 1) is an arbitrary matrix-valued function such that 



3,0). z, ffl 



= z. 

f=0 



Statistics of energy partitions in ai'bitrary dimension 



7 



In all the cases below, 3£ will be well defined for all the pairs (Z, Z) with the possible ex- 
ception of matrices Z lying in a set of positive codimension. For such Z, one will always be 
able to define the resulting quantities we will be interested in by continuity. 



3 Hyperangular momenta 

Instantaneous phase-space invariants dependent on Z are exemplified by the so-called hyper- 
angular momenta of a system of particles (with respect to the origin), namely, the physical 
angular momentum J, the kinematic angular momentum K dual to 7, the grand angular 
momentum A, and the singular angular momentum L (8 l l9l llOlfTT]|13l . These non-negative 
quantities are defined by the formulas 

/2= 52 Jij=Mf^{ZiaZja-ZjaZia), (6) 

\<^i<i!id a=\ 

d 

= 52 ^a^' = Mj^{ZiaZip -ZipZia), (7) 

A^=m' 52 (Z,„Z,;j-Z,,jZ,„)2, (8) 

i<j or i=j&.a<p 

L^= 52 Ll„ La,=M{^^i,-^,^a) (9) 
l^(T<T^m 

(^1 , ^2, • • • , "^m being the singular values of Z and their derivatives E,a being defined in Re- 
mark[TJ. Of course, the formula (|6) coincides with (|2) if Z is the position matrix in a Carte- 
sian coordinate frame and n = A^. In our previous papers [8 9 10 13], we denoted the singu- 
lar angular momentum by . 

The correct definition Q of the kinematic angular momentum K was first given in our 
article ||9l in the particular case d = 2> although the words "kinematic angular momentum" 
themselves were used earlier fT'Sl. The grand angular momentum A was first introduced by 
Smith |2|, also in the particular dimension d = !> (and mainly for systems of three particles). 
The singular angular momentum L was first defined in the paper |[8], again in the particu- 
lar case d = 3 only. For an arbitrary dimension d of the physical space, the hyperangular 
momenta J, K, A, and L were introduced in our article 0101 . 

Theorem 1 All the four hyperangular momenta J, K, A, and L are instantaneous phase- 
space invariants. 



= E K'rjJ)-rjJ)r£)l (lO) 



Theorem 2 The hyperangular momenta J, K, and A can be alternatively computed as 

'^ap'^ap '^ap'^Pa 

a.p=\ " 

(1) ^' (2) ^ ■ (3) . . 

/=1 /=1 

-M^i (11) 
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A7 — ^ ZiaZja, A- J — Y^ZiaZja, A^j — ^ ZiaZja, 
a=l a=l a=l 

A^=M^(\\Zf\\Zf-{Z,Zfy (12) 
There hold the inequalities +1? ^ and ^ A^. 

Remark 2 One of the important consequences of Theorem |2] is as follows. The formulas 
©-([oil imply that the numbers of components of the hyperangular momenta J, K, A, and L 
are equal respectively to d{d — l)/2, n{n — l)/2, dnidn — l)/2, and m(m — 1) /2. However, 
the alternative formulas dllt and (I12t show that K and A can be calculated (like J) using no 
greater than operations where the constant Q depends on the dimension d only. 



4 Energy partitions 

The instantaneous phase-space invariant 

r = -|izf (13) 

is called the total kinetic energy (of the system of particles described by the position matrix 
Z and its time derivative Z). Of course, the formula (I13t coincides with ([TJ if Z is the time 
derivative of the position matrix in a Cartesian coordinate frame and n=N. 

The main subject of this article is statistical properties of the terms of various partitions 
of the total kinetic energy T (almost all of these partitions were introduced in our paper (9) 
for = 3 and in the subsequent paper ifTOl for an arbitrary dimension d of the physical 
space). In order to define the partitions in question, one has to make oneself more acquainted 
with the action ^ of the group 0{d) x 0(n) ifTOl . 

Given ad x n matrix Z, consider the manifolds of all the matrices of the form 

{RZ\ReSO{d)}, {Ze| eeSO(n)}, {RZQ \ R e SO{d), Q e SO{n)}. 
Denote respectively the tangent spaces to these manifolds at point Z by 

nieft(z) = {z+3?z 1 3?eso(d)}, 

nright(Z) = {Z + ZQ I Qeso(n)}, 

n{Z) = {Z + JiZ + ZQ\5le30{d), Q6so(w)} 

[recall that 3o{d) and so{n) are the spaces of skew-symmetric d xd and nx n matrices, re- 
spectively]. The multidimensional planes nieft(Z), nright(Z), and n(Z) are affine subspaces 
of the space W^^" of d x n matrices. 

Theorem 3 The intersection of the spaces J7ieft(Z) and /7right(Z) consists of the only matrix 
Z: 

Aeft (z) n nrigh, (z) = {z} ^ dim n^t (z) + dim nrfght (z) = dim n (z) 

^ n (z) - z = [Aeft (z) - z] ffi [Hright (z) - z] 

if and only if all the positive singular values of the matrix Z are pairwise distinct. 
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Given now adxn matrix Z, denote by Z'*^^', Z"^''', and Z^°' the orthogonal projections of 
Z onto the spaces JTieft(Z), J7right(Z), and TI{Z), respectively, in the sense of the Frobenius 
inner product ([5). Let also Z' = Z.—Z}°^ denote the component of Z orthogonal to n{Z). 

Introduce the following non-negative quantities: 





A2 




Ta 




called the grand angular energy, 




Mp2 
2 


called the hyperradial energy, 


y^rot 




called the rotational energy. 


T' 




called the inertial energy. 




2Mp2 


called the shape energy. 



Ta, Tp, and T^ being defined for p > only (concerning p, see RemarkO. 
Theorem 4 There hold the identities and inequalities 

T = TA+Tp, (14) 

r = r°' + r^ (i5) 

r"' s^Ta^^T, Tpi^T' T, Ta - T'°^ = t' -Tp = T^^ mm{TA,T'), 

Definition 3 The equality T = TA+Tp\& called the Smith decomposition of the total kinetic 
energy T . The equality T = 7'°' + T' is called the orthogonal decomposition of the total 
kinetic energy T . 

The Smith decomposition lll4t was deduced for the first time by Smith 1 2 1 in the partic- 
ular case N = d = 3. 

The identity ( 115b is in fact just the Pythagorean theorem: it follows immediately from 
the relations Z = Z'^^+Z' and (Z'°',Z^) = 0. The identity GUi is an immediate consequence 
of the equality ( 112b . Indeed, since p = ||Z||, one has pp = (Z,Z). Now 

2T _ ^ 2_ \\Z\\HZ\\'-{Z,Z)^ , {Z,Zf ^ ,_2Ta 2Tp 

M " " p2 ^ p2 M2p2 +^ M ^ M ■ 

Introduce now the non-negative quantities 

jext ^ ^ I i^left 1 1 2 called the external energy, 

r'"' = ^ I |Z"sht 1 1 2 called the internal energy. 



2 

72 

7} = ^ called the outer angular energy, 

Tx = J called the inner angular energy 
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as well as the quantities 

T^ies ^ ji-ot _ T^ext _ T^int ^^jj^j jj^g residual energy, 

= 7"°' — Tj — Tk called the angular coupling energy, 

Tj, Tic, and being defined for p > only. The energies r"'" and Tac can be negative. 
Definition 4 The equality 

j^rot j^ext ^ j^int _j_ j^ies (16) 

is called the projective partition of the rotational energy T™'. The equality 

r°' = Tj + TK + T.,, (17) 
is called the hyperspherical partition of the rotational energy 7"'"°'. 
Tiieorem 5 There hold the inequalities 

Tj + T^^Ta, Tk + T^^Ta, 

Tj ^ 7"="' s; r™' , Tk s; r'"' ss r°' , r^' «; Jac ■ 

If d ^2 or n = \ then Tj = r""'. Recall that the case n = 1 corresponds to the situations 
where either there is only one particle {N = 1) or the system consists of two particles {N = 2) 
whose center-of-mass is fixed and coincides with the origin. Ifd = \ or n ^2 then Tk = T'"'. 
Recall that the case n = 2 corresponds to the situations where either there are two particles 
(N = 2) or the system consists of three particles {N = 3) whose center-of-mass is fixed and 
coincides with the origin. 

Finally, suppose that all the positive singular values of the matrix Z are pairwise distinct 
(besides these positive singular values, the matrix Z is allowed to possess zero singular value 
of arbitrary multiplicity). Then, according to Theorem |3] the component Z™' £ [n{Z) — Z] 
of Z can be uniquely decomposed as 

z"' = z™'+z'", z™' 6 [nieft(z) -z], z- e [nn^z) -z]. 

Introduce now the non-negative quantities 

j^oui _ ^ |j20ut|j2 called the tangent (or singular) external energy, 

E'" = — ||Z'"||^ called the tangent (or singular) internal energy 

as well as the quantity 

^coupl _ jrot _ ^oiit _ ^in called the tangent (or singular, or Coriolis) coupling energy 

which can be negative. In our previous papers l|9l|T0l. the quantities E°"\ and £''=°"Pi 
were called respectively the outer term, the inner term, and the coupling. The energy £''^°"pi 
was called the Coriolis coupling term in ||5J (see also |i6J). 

Definition 5 The equality 

r™'=£™' + £'"+£"™P' (18) 
is called the singular value expansion (or just singular expansion) of the rotational energy 

j^rot 
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The singular value expansion dlSt is defined only if all the positive singular values of 
the matrix Z are pairwise distinct. However, this condition is met for generic matrices Z £ 
j^dxn |jQj Some partitions of the total kinetic energy T essentially equivalent to the equality 
T = ii™' + £■'" +£''^°"pi + t' have been known long ago (see e.g. the papers |[Tj|3])- 

Theorem 6 Let Z = DYX* be a signed SVD ^ of a d x n matrix Z. Suppose that all the 
three derivatives D, T, and X are well defined (see Remark\J}. Then 

Deso(rf), Xeso{n), Z = bTX* +DtX* +DTX\ 
Z'°* = brx* + DYX* , Z' = DYX* , 

{z+DYx*) e n,eft(z), {z+DYX*) 6 nright(z). 

In particular, if all the positive singular values of the matrix Z are pairwise distinct then 

Z°"'=bYX*, 2'"= DYX*. 

Theorem 7 All the 14 energy quantities Ta, Tp, T™', T', T^, T'^"', T'"', Tj, Tk, T^c, 
i?""', £■'", and £'^°"Pl introduced in this section are instantaneous phase-space invariants. 

For a fixed total kinetic energy T, the energies J^, Tp, J™', T' , T^, J*^"', 7'"', Tj, Tk, 
T""^^, and [the terms of the Smith decomposition (114b . the orthogonal decomposition I llSt , 
the projective partition ( I16t . and the hyperspherical partition jnH are bounded according to 
Theorems |4] and |5] 

q^Ta^t, os^Tpi^T, ^ 7"'°' ^ r, ^ r' ^ r, o ^ ^ r, 
^ r""' ^ r°' , ^ r'"' ^ r"' , | r''' \ ^ , 
Osj7>^r"', os;7xr°', |7;c| r"'. 

On the other hand, the terms ^ 0, > 0, and £'=°"Pi T"' of the singular value 
expansion ( 118b are no?: S""', and — £<^°"p' can be arbitrarily large for any fixed value 
r > HOlliri . That is why we call the equality dlSb an expansion rather than a partition and 
use the letter E for its terms rather than the letter T. 

Remark 3 In the physically most important case where Z is the position matrix of a sys- 
tem of A' particles (a cluster) whose center-of-mass is fixed at the origin, the meaning of 
all the terms Ta, Tp, r°\ T' , T^, T^''', T"', 7>, Tk, T'^', T^^, E""', and £™"p' is dis- 
cussed in detail in our previous papers ll9] |101fT3l (where also the names of these terms and 
the words "orthogonal decomposition", "projective partition", "hyperspherical partition", 
"singular value expansion" are justified). To be brief, the hyperradial energy Tp corresponds 
to the contribution (to the total kinetic energy T) of changes in the size of the cluster as a 
whole and the shape energy T^ , to the contribution of changes in the shape of the cluster. 
The inertial energy T' = Tp + T^ describes the contribution of all the changes in the singular 
values i§i,(§2,...,i§(B where CO = mm{d,N — 1). The external energy T'^"' is the contribu- 
tion (to the total kinetic energy T) of rotations of the cluster in M.'^ (about the origin) as a 
whole. These rotations leave the symmetric N xN matrix Z*Z invariant. The internal energy 
r'"' is the contribution of so-called ''kinematic rotations" Il9l ll0llllll35l of the cluster in the 
"kinematic space" R^^' as a whole. "Kinematic rotations" are the cluster rearrangements 
that leave the symmetric d xd matrix ZZ* invariant. The rotational energy 7'°' corresponds 
to the contribution of rotations of both types, physical ones and "kinematic" ones, and the 
residual energy T'"*^^, to the coupling between these two types of rotations. Finally, the grand 
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angular energy TX = 7"'°' + describes the joint contribution of conventional rotations of 
the cluster, "kinematic rotations", and "rotations" of the vector (i^i , (§2, ■ • • , i^o) of the singu- 
lar values in R'". The physical meaning of the terms of the hyperspherical partition ( 117b and 
the singular value expansion ( 118b is much less clear. 



5 Refinement of the singular value expansion 



The singular value expansion dlSb of the rotational energy 7""* can be partitioned further as 
follows 1 11]. Suppose that among the singular values £,1,^2,- --i^m of the matrix Z, there 
are A: ^ m positive numbers, and these numbers are pairwise distinct (which ensures the 
existence of the singular value expansion): 

<§1 > <§2 > ^k+\ = • • • = <§m = 0. 

Then the symmetric d xd matrix ZZ* and the symmetric n x n matrix Z*Z have k positive 
pairwise distinct eigenvalues £,f,£,2,...,£,l, whereas the remaining eigenvalues (if any) of 
each of these matrices are equal to zero. We will denote the standard inner product in Mf' 
and M" by x ■ y and the Euclidean norm (x • x)'/^ of a vector x by |x| (in fact, we already 
used the notation |x| in Sections [T] and |2]l. 

Denote by U\,U2, ■ ■ ■ ,iik the unit eigenvectors of the matrix ZZ* corresponding to the 
eigenvalues £,f,£,2, • ■ • , respectively. These vectors are determined unambiguously up to 
multiplication by —1. Let Uq be the zero subspace (of dimension d — k) of the matrix ZZ*, 
i.e., [/q = {u 6 M'' I ZZ*u = 0}. Denote by a,- the orthogonal projection of the derivative li/ 
of the vector u, (see Remark[Tli onto Uq in the sense of the standard inner product x • y in R'' 
(1 ^k). Introduce the non-negative quantities 

M * 

= ^ ^ "3/ 52 ' called the unbounded component of the outer term, 

2 ,=1 >=i 

M 

= -z-'^^ll^if' called the bounded component of the outer term. 

2 1=1 



Similarly, denote by vi , V2, . . . the unit eigenvectors of the matrix Z*Z corresponding to 
the eigenvalues ^f,£,2,...,^^, respectively. These vectors are again determined unambigu- 
ously up to multiplication by — 1 . Let Vq be the zero subspace (of dimension n — k) of the 
matrix Z*Z, i.e., Vq = {v £ R" | Z*Zv = 0}. Denote by ha the orthogonal projection of the 
derivative \a of the vector \a onto Vq in the sense of the standard inner product x ■ y in M" 
(1 ^ a ^ k). Introduce the non-negative quantities 

£|" = ^ 52 ''a H (^|6 '^a)^ called the unbounded component of the inner term, 

a=\ ;6=1 




called the bounded component of the inner term. 



It is clear that the energies ii™' and iij do not change if one multiplies some of the vectors 
Ui,U2, . . . ,U4 by —1 and the energies £[" and is™ do not change if one multiplies some of 
the vectors vi , V2 , . . . , by — 1 . 
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Remark 4 It is obvious that u, • Uj = — uy • li, and u, • li, = for any 1 ^ i,j ^ k. Conse- 
quently, 

i:^2^(u,.u,-)^ = i f:(42 + ^J)(u,.u,)2= J: & + ^f){nr»j)\ 

1=1 j=i u=i K/<ja 

so that 

Er = ^ E {^l + ^j){nrnj)\ (19) 

and analogously 

£r = y E (^a + ^p')(Va-V/if. (20) 

For an arbitrary dimension d of the physical space, the terms Z?™' [in the form (I19H. isj" 
[in the form ( i20b 1, fi™'' ^2° were introduced in our previous paper Hill . However, in the 
particular case d — 3, the energies and were defined in the literature earlier (using 
other names). For instance, the formula (33) of the article [5| expresses the total kinetic 
energy T of a cluster of A' ^ 5 particles as the sum of five components. In our notation, 
those components are ^^"p', T', and Ef. 

Remark 5 Let Z = DTX* be a (signed) SVD Q of the matrix Z with 

Then u, can be chosen to be the i'th column D, of the matrix D for each I ^ i ^ k. Indeed, 
ZZ* = DTY*D* . If e/ denotes the /th unit coordinate vector of R'', then 

ZZ*Di = DTT*D*Di = DTT*ti = E,fDti = E,fDi. 

Similarly, \a can be chosen to be the ath column Xa of the matrix X for each \ ^ a ^k. 
Indeed, Z*Z = XT*TX*. If ta denotes the ath unit coordinate vector of R", then 

Z ZXa = XT TX Xa = XT Tta — "ra-^Cct = ^^a^a ■ 

Theorem 8 All the four additional energy quantities Zi°"', E^^^, and E™ are instanta- 
neous phase-space invariants, and there hold the decompositions 

£out^£out_^£,out^ £'" = £;"+£:|". (21) 

Moreover, the inequalities 

^out ^ j^rot ^in ^ j^rot 

are valid. 

In our previous paper Hill , this theorem was just announced, so we will prove it here. 

Proof Examine the unbounded and bounded components E™ and E™ of the inner term. 
Consider a transformation Z' = RZQ* , Z' = RZQ* with/? 6 0{d), Q e 0(«). Then (Z')*Z' = 
QZ*ZQ* and 

Vi = QVq, 

v'„ = ev„, v'„ = ev„, b'„ = eba (i^ 
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Consequently, the quantities Zi|" and Zij" remain unchanged after such a transformation. If 
one augments both the matrices Z and Z by the (« + l)th column equal to zero: Z' = (Z 0), 
Z' = (ZO),then 



(n+l)x(n+l) 



3 (Z')*Z' : 



Z*Z 




Thus, 



Vq' = {(vw) e R"+' I V e Vo, w 6 R}, 
v'„ = KO), v; = (v„0), b'„ = (baO) {l^a^k), 

and the quantities £[" and is™ again remain unchanged. Thus, these quantities are instanta- 
neous phase-space invariants. 

The relations £j" -|- Zij" = and Ef ^ 7™' will be verified only in the case where the 
pair (Z, Z) admits a signed SVD Z = DYX* of the matrix Z with T as in Remark |5] and 
with the well defined derivatives D, T , and X (see Remark[T]l. The opposite case is highly 
degenerate and of no practical importance; in fact, the equality Zsj" + Zi™ = Z?'" and the 
inequality E™ ?C J™' will follow for the exceptional pairs (Z,Z) by continuity. 

According to Remark|5] one can set \a = Xa and therefore \a =Xa ^ oc ^ k), where 
Xa denotes the ath column of the matrix X and Xa denotes the ath column of the matrix X. 
According to Theorem |6] Z'" = DTX* and Z"' = DTX* +DTX*. Since Ej", Ef, and 
T"' are instantaneous phase-space invariants, we can assume without loss of generality that 
D is the identity d x d matrix and X is the identity nx n matrix. Then 

z'" = rx* , z'°' = or + rx* , 

Va = Ca is the ath unit coordinate vector of R", 

rz = diag(^2,^|,...,^2^ 0,0,. ..,0), 

n-k 

and Vf) is the subspace of all the vectors v of the form 

v= (0,0,...,0,wi,W2,...,w„_,t). 
k 

Thus, 

V|6 -Vce =Xpa, ba = {0,0, . . . ,0 ,Xk+iM,Xk+2.a, ■ ■ ■ ,Xna) 

k 

(1 ^ 0!,/3 ^ fc). Consequently, 

= y L L ^/3a> Ef = —Y,^a L ^j3a' ^^^^ 
a=l 13=1 a=l j3=it+l 



and 



^ a=l |6=1 ^ ^ 

Moreover, the last w — A: columns of the dxn matrix t)T vanish. Therefore, 
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is no less than M/2 multiplied by the sum of the squares of all the entries of the last n — k 
columns of the x n matrix TX* , i.e., is no less than Ef [see ( I22b 1. 

The unbounded and bounded components Z?™' and £2"' of the outer term can be treated 
the same way. Note only that if one augments both the matrices Z and Z by the (n+ l)th 
column equal to zero: Z' = (Z 0), Z' = (Z 0), then Z'(Z')* = ZZ* , Ul^ = Uq, and uj = u,-, 
uj = u,-, aj = a,- for each I i^i ^k. □ 

However, it is rather difficult to reveal the physical meaning of the quantities Zi™', Z?™', 
Ef, and Ef. Of course, ii™' and E'" can be arbitrarily large for any fixed value T > (like 
£™' and 



6 Equal masses and random masses 

To study the statistical properties of various components (defined in Sections |4] and [Sjl of 
the total kinetic energy T of systems of classical particles in R"', one has to describe pre- 
cisely the sampling procedure for the coordinates and velocities of the particles. We will 
consider systems of N ^2 particles with the center-of-mass at the origin. As in our previous 
paper [11], two situations will be dealt with: particles with equal masses and particles with 
random masses. 

Let Wa and Wa (1 ^ 0! ^ A') be IN independent random vectors, each being uniformly 
distributed in the unit ball in centered at the origin. Now set 

(the subscript "cm" is for "center-of-mass"). The masses ma of the particles in the situation 
of equal masses are equal to 2/N: 



N 

mi = /«2 = • ■ • = = — , M= Y^mc 



and one computes 

ga=Wa-Wcm, ga=Wa-Wcm, l^a^N. 

In the situation of random masses, one first chooses the masses ma of the particles according 
to the formula 

( N N 

m« = 2r]„ E , l^ai^N, M = E '«« = 2, 

V/S=i / «=i 

where r]i , r]2, • . • , 7]^ are independent random variables uniformly distributed between and 
1. The vectors ga and ga are then computed as 

ga =ma'''^(Wa-Wcm), ga = '«a'''^(Wa - Wcm) , l^a^N. 

In both the situations 

E '«« ga = E "'a gcf = 0. (23) 

a=l a=l 
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Finally, the mass-scaled radii vectors and their time derivatives are calculated in both 
the situations as 

qa = Ciga, qa = Giga, 1 ^ « ^ A^, (24) 
where the positive factors Ci and 62 are determined from the condition 

Ilqa|'= E |qa|' = l- (25) 

a=l a=l 



In both the situations 



E ™a ^qa = E ™« = (26) 



according to i23i . so that the d x N position matrix Z with columns qi , q2 , • • • , qA' and its 
time derivative Z with columns qi , q2, . . . , qM do correspond to a system of particles with 
the masses m\,m2,. .. ,111^ and with the center-of-mass at the origin (see Example [T). More- 
over, the normalization ( 125 1 is equivalent to that ||Z|| = ||Z|| = 1. Together with M = 2 this 
gives p = T = I. 

Remark 6 A random point w uniformly distributed in the unit ball in R"' centered at the 
origin can be generated as w = x'/'^s where x is a random variable uniformly distributed 
between and 1 and s is a random point uniformly distributed on the unit sphere S'^^ ' C Mf' 
centered at the origin (h and s being independent). There are many methods to choose s for 
d ^2, see e.g. the papers |41 42 43 ,44,45 1 and references therein. The most "elegant" (and 
probably the best known but not the fastest) algorithm is the so-called Muller or Brown- 
Muller procedure f42| which runs as follows. Let XiyX2,- --tXcI be independent random 
variables normally distributed with zero mean and the same standard deviation g > 0. Then 
the point 

-(y'f ¥)■ '^-p <»^»»- 

is uniformly distributed on S''^'. Indeed, the joint probability density function 



t.2 



oiX\,X2,---,Xd depends on 5 only. 
Of course, for <i = 2 one can just set 



s = (cos(p, sin(p), 



where 9 is uniformly distributed between and 2n. For d = 2> the standard spherical coor- 
dinates in M? lead to the choice 111] 

S= ((l-/;2)l/2cos(p, (1-/72)1/2 sin 

where h and 9 are independent random variables uniformly distributed in the intervals — 1 ?C 
h ^ 1 and ^ 9 ^ 2;r. 

The following conjecture is the main statement of this paper. 
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Conjecture 1 Introduce the notation 

v=N-\, co = mm{d,v). 

For any d ^ I and N ^ 2, the mathematical expectations E of the bounded terms 7X, Tp, 
r°\ T', T^, 7"="', r'"', Tj, Tk, T^c, ^r'- and of the Smith decomposition 

the orthogonal decomposition ( 115b . the projective partition l ll6b , the hyperspherical parti- 
tion ( 117b , and the singular value expansion dlSb , ( 121b in the situation of equal masses are 
given by the formulas 

¥.Ta = \-—, Er„ = — . Er°' = l-— . ¥.T' = —, (27) 
dv ^ dv' dv' dv 

Er,=^, (28) 

^^^^^^cojld-co-l) ^^,.^c.(2v-a>-l)^ 
2rfv 2dv ' 

d-l v-1 J + v + a)-2 

E7>=-— , Er^=-— , E7;c = l , (30) 

dv dv dv 

EET=1--, e4" = 1--. (31) 
d V 

Since we use the normalization 7 = 1 , the formulas (I27b-(l31b give in fact the mathemat- 
ical expectations of the ratios T^/T, Tp/T, r°*/T, T'/T, T^/T, T'^'^T, P'^^/T, r'^'/T, 
Tj/T, Tk/T, T^c/T, E^^^/T, and E'^"/T. These formulas reflect vividly the duaUty of the 
physical space R'' and the "kinematic space" M*' (see Example[T]l. The expression for each 
of the quantities ET^, ETp, Er°\ ET', ET^, E^"=^ and ET^^ is symmetric with respect to 
d and v. If one interchanges d and v, the expressions for ET"™', ETj, and Eii™' turn into 
those for EJ'"', ETk, and EZsj", respectively, and vice versa. 

Note that all the four formulas ([27}, (I2S}, and ([30} for the mean values ET^, ET^, E7>, 
and ETis: of the energy terms H corresponding to the momenta A,L, J, and K have the form 

7- 1 

EH=^, 
dv 

where, roughly speaking, y is the dimension of the space where the "rotation" responsible 
for the momentum in question takes place. This space is the matrix space W'^^ for A, the 
space { ((§1 , ^2, • . • , ^(o) } of the collections of singular values for L, the physical space M'' 
for /, and the "kinematic space" W for K. 

Example 3 Consider a system of two particles (A' = 2) of arbitrary masses m\ and m2 with 
the center-of-mass at the origin of W' . For = 2 it is not hard to obtain IIOIIIII that 

m\ mi , ,T 
2M ' ' ' 

_ yint . j^res _ ^in ^in ^in ^^^^ ^coupl q 

Tp = t' = Tcos^d, 

= r°' = r"' = Tj = E°'" = = rsin^ e, 

where r is the vector connecting the particles and Q is the angle between r and r (cf. Sec- 
tion [T]l. It is well known 111 J that 

Ecos^ B = \/d, Esin^e = (t/- 1)/J (32) 
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for the uniform distribution of the directions of the vectors r and r in R"'. One thus concludes 
that for = 2 and r = 1 

Er^ = EJ'"' = Er'''=' = ¥.Tk = ET^ac = E£'" = E£j" = E^" = ESj"" = E£™"p' = 0, 

Erp=Er' = i/j, (33) 

KTa = E7""°' = Er'^"' = E7> = E£°"' = EEf"' = {d-\)/d (34) 

independently of the masses m\ and mi- This agrees with the formulas (I27t -(l31t for O) = 
V = 1. Thus, for A' = 2 Conjecture [T]is correct for any d (and, moreover, the equality of the 
masses of the particles is irrelevant for A' = 2). 

The formula (|27} for ET"' and EJ^ as well as the formula (|29} for E7''='", EJ'"', and 
Er"^'* were proposed in our previous paper 0111 . Moreover, the paper 111] contained also 
the equalities (I27b-( l31b for all the terms in the particular dimension d = 3. The formulas 
(I27b-( l31b for d = 3 were confirmed in the paper 0111 for 3 A' ^ 100 by extensive numerical 
simulation. Note that for d = 3, the equality EE™' = I — co/d reduces to 

2/3 for N = 2, 
1 /3 for N = 3, 
for A? ^ 4. 

The equality Eii™' = for d = 3 and N ^4 (independently of the distribution of the masses) 
is obvious [see i35l below]. The numerical experiments of the paper 1111 for f™' were of 
course carried out for A' = 3 only. 

Conjecture 2 The equality ET"^''*' = holds for any "reasonable" distribution of the masses 
of the particles (not only in the situation of equal masses). Moreover, T'''^ > with proba- 
bility 1 /2. 

Conjecture |2] was confirmed in the numerical simulation of our previous paper Bill for 
d = 3 and 3 ^ A^ ^ 100 in the situations of equal and random masses. Some "theoretical" 
support for this conjecture was also presented in fll]. 

7 A "physical" proof of Conjecture[T] 

Consider a system of N ^ 2 classical particles in R'' with masses mi,m2, . . . ,mff and with 
the center-of-mass at the origin. Such a system can be described either by the if x A^ position 
matrix Z (in a Cartesian coordinate frame) and its time derivative Z with the columns subject 
to the identities ( I26l l or by d x (A' — 1 ) matrices Z and Z with arbitrary columns (see Ex- 
ample O. We will follow the second approach and regard Z and Z as matrices of Frobenius 
norm 1 [cf. ( I25b 1 with d rows and N — I = v columns. Recall that M = 2 according to the 
procedures of Section |6] 

The following lemma is well known llll and almost obvious. 

Lemma 1 Let a random point x be uniformly distributed on the unit sphere Ef^''^ ' C R'''^'' 
centered at the origin (p, q £ N), and let be an arbitrary affine plane in R^+* of dimension 
p. Denote by Pr|,x the square of the length of the orthogonal projection ofx onto <i>. Then 
the mathematical expectation EPr^ x o/Pr^ x is equal to p/{p + q). 



E£|"' = 
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Proof Let {xi,X2, . . . ,Xp+q) be a Cartesian coordinate frame in R''+'?. Without loss of gen- 
erality, one can assume that the plane is given by the equations x^+i = = • • • = 

Xp+g = 0. If X = (xi,X2, . . . ,Xp+g) e gP+i-^ then Pr^x = x]+xl H \-xf,. But xj+xj + 

\-x^_^_g = 1 and E;C[ = Ex2 = • • • = ]E;c^,_|_^ by symmetry reasons. Therefore, E,x^ = Kx2 = 

■ • • = Exj+g = 1 / (p + <?) and E Pri X = /?/ (p + 9) . □ 

By the way, the formula i l32t for Ecos^ is a particular case of Lemma [T] for p = I, 
q = d — I. 

The following reformulation of Lerrmia[T]is especially useful: if a random point x = 
{x\,X2, ■ ■ ■ ,Xp+q) is uniformly distributed on the unit sphere S''^''^^ centered at the origin, 

then E(xjj H i"^? J — p/iP + l) fo'' ™y fixed set of indices 1 ^ ii < I2 < • • • < ip ^ 

p + q. 

All the CO = min(J, v) singular values of a generic d x v matrix Z are positive and 
pairwise distinct, whence (see LIO]) 

„ a(2d-co-l) „ co(2v-co-l) 
dimnieft(Z) = ^ ^, dimnright(Z) = ^ 

dim/7(Z) = dimnieft(Z) -|-dimnrigh,(Z) = dv - CO 

[note that co{d + v — o) = comax{d, v) = dv]. Let Z be a random matrix uniformly dis- 
tributed on the unit sphere S''*'^^ c R"'^*' = K"'^ centered at the origin [observe that after 
the identification R''^^ = M''^, the Frobenius irmer product (|5} becomes the standard inner 
product]. Since M = 2, 

j^ext ||2^1eft||2 j'int ij^rightjp j^rot ||^rot||2 

On the other hand, Z''^*', Z"sht^ gj^^^ j^ox orthogonal projections of Z onto the spaces 

nieft(Z), n,-ight(Z), and n(Z), respectively, in the sense of the Frobenius inner product. 
According to Lemma [T] 

^^,,^co(2^-a)-l) ^^.„.^a>(2v-a)-l) ^7- = !-^, 
2dv 2dv dv 

and 

Er^ = 1 - EJ"' = , EJ''"^ = Er"' - Er""' - Er'"' = o. 

dv 

We have verified the formulas ([27} and (|29li for EJ™', EJ', EJ'^'", EJ"", and £7''"'='* follow- 
ing the reasoning in our previous paper [ill j . 

Since p = ||Z||, one has p = ||Z||^' (Z,Z) = ||Z||cosiJ, where •& is the angle between 
Z and Z in the sense of the Frobenius inner product. Let us treat Z and Z as independent 
random matrices uniformly distributed on the unit sphere in R"'^*' centered at the origin. For 
Tp = Mp2/2 = p^ = |lZ|p cos- ^ = cos^ ■& we therefore have 

Erp = -^ 

^ dv 

[see ( 1321 : the equality ||Z|| = 1 is in fact irrelevant here], and 

EJa = l-EJo = 1-— , E7> =Er^-Er„ = 

dv ^ ^ dv 

We have verified the formulas ([27) and ^ for EJa, ET^p, and EJ^. 
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To calculate ETy and E7)f , it is expedient to use the expressions (do} and ([IB for and 
K^, respectively. Since M = 2 and p = 1, 

a,;8=l 

where 

(1) (2) . (3) . . 

/=1 /=1 /=1 

Let again Z and Z be independent random matrices uniformly distributed on the unit sphere 
in R''^^ centered at the origin. Then 

whenever a ^ ji. Indeed, if one changes the signs of all the entries Z/a, I ^d, the sums 
r^^p and r^^^ would change their signs while the sums and F^'^' would remain the 
same. Consequently, 

ETy = E f {r^'Jr^'d - frj^f I = e f [ri^Vj^' sin^ 0„ 

= E EFii'Erj^^Esin^ 6. = = 

where Set denotes the angle between the ath column of Z and the ath column of Z in R'' 
(1 s£ a s; V). Here we have used the facts that for each cu, the three random variables r^a 
(the square of the length of the ath column of Z in R''), F^a (the square of the length of the 
ath column of Z in R''), and sin^ Ga are independent, and according to Lemma [T] 

Erii)=Eri^' = A = l, Esin2e„ = ^ 

dv V d 

[see (l32li1. Analogously, 

V- 1 

ETk=-—, 
dv 

and 

ETLc = Er°' -ETj-ETk = 1- ^ + ^ + ^~^ . 

dv 

We have verified the formulas (I30t for E7), EJ/^ , and EJ^c- 

Finally, consider Zij"' and Zs™. Generically, in the notation of Section [S] 

k = d, Uq = 0, E^"' = 0, £™' = £™' for di^v, 
k=v, Vb = 0, 4" = 0, for v J. 



(35) 



Let V > J and find EZ?"- Assume all the d singular values ^1,^2, ■ ■ ■ of the matrix Z to be 
positive and pairwise distinct (a generic setup). Then in the SVD Z = DYX* of the matrix 
Z, all the three derivatives D, Y, and X are well defined, and 

Z = DYX* + DYX* + DYX* (36) 
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according to Theorem |6] Since Z and Z are independent and Zij" is an instantaneous phase- 
space invariant (Theorem [8), one may suppose that D is the identity d •>< d matrix, X is the 
identity v x v matrix, while Z is still a random matrix uniformly distributed on the unit 
sphere in M''^^ centered at the origin. The equality (I36t takes the form 

Z = t)T + t + TX*. (37) 

The last v — d columns of the matrices DT and T vanish. Consequently, the equality ( I37l i 
implies 

for any d+l ^ a and I ^i^d. Taking into account that M = 2, the second formula i22l 
[valid for identity matrices D and X] becomes 

d V d V 

1=1 a=d+l i=l a=d+l 

According to Lemma[T] 

„ in d(v — d) d 

EEf = ^^ = l--. (38) 

Now observe that the formula (I38t for v > J and the equality Eiij" = for v d can be 
combined into the unified formula 

e4" = 1 - - 

V 



valid for any d and v. Analogously, 



E£2°"' = 1 - 



d 

for any t/ and v. We have verified the formulas dlD for E£f"' and Efif- □ 

Unfortunately, all these arguments are not mathematically rigorous because even in the 
situation of equal masses, the procedure of generating the matrices Z and Z described in Sec- 
tion|6]does not give (after the passage from W'^'^ to M''^('^^') = M'''^^) matrices uniformly 
distributed on the unit sphere in Mf^^^ centered at the origin. This can be easily shown even 
in the simplest case d = I, N = 3. For these values of d and A', the procedure of Section|6] 
for equal masses (in the case of Z for definiteness) runs as follows. 

Let w\ = a, W2 = b, = c he three independent random variables uniformly distributed 
between —1 and 1 (in the "unit segment" in R). One computes the numbers 

2a — b — c 2b — a — c 2c — a — b 
g] = ^ , 82 = ^ , g3 = ^ 



and the vector 

G{2a — b — c, 2b — a — c, 2c — a — b) 



of length 1 [in the notation of ll24t . C = C2/3 > 0]. Now one has to choose an arbitrary 
matrix 

/ Gil G12 Gi3 \ 

G = G21 G22 G23 e 0(3) 

1 3-1/2 3-1/2 3-1/2 / 
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(see Example O and calculate the vector 

e{2a-b-c, 2b-a-c, 2c-a-b)Q* . (39) 

The first two components of the vector (I39t are 

Zi=e[{2a-b-c)Qn + i2b-a-c)Qi2 + {2c-a-b)Qi3], 
Z2 = e [(2a -b-c)Q2i+{2b-a- c)Q22 + {2c -a- b)Q2i] 

{Z\ +Z| = 1), the third component is zero. Then (Z\,Z2) is the 1x2 matrix Z one deals 
with. 

Do there exist fixed numbers Qw, Q12, Qu, Q21, Q22, Qn such that the point (Zi,Z2) 
is uniformly distributed on the unit circle §^ C centered at the origin? The answer to this 
question is negative. Indeed, let A > 0. The inequalities 

0^Zi^AZ2, (40) 

that is 

sS (2a - ^7 - c)en + (2& - a - c)ei2 + (2c - a - fo)ei3 
s: A [(2a - ^> - c)e2i + (2fc - a - c)e22 + (2c - a - b)Q2-i\ , 
-1 ^ a ^ 1, 
-1 ^^-^ 1, 
-1 SjcsC 1, 

determine a certain polyhedron in the Euclidean space M.^ with coordinates a, b, c. The coor- 
dinates of the vertices of this polyhedron and its volume V are piecewise rational functions 
of A, Qn, Q\2, 2i3, 22b Qii, G23 with integer coefficients (see Remark |7] below). Conse- 
quently, the probability V/8 of the inequalities l l40b is also a piecewise rational function of 
■^i Gill Gi2i Gi3> G21, G22> G23 with integer coefficients [the denominator 8 in V/8 being 
the volume of the cube max(|a|, |fe|,|c|) ^ 1]. On the other hand, if the point (Zi,Z2) were 
uniformly distributed on S', then the probability of ( 1401 ) would be equal to ^ arctanA. 

Remark 7 Consider a collection of DC affine hyperplanes 

yi 

7=1 

in R-'^ = {(xi,X2, . . . jXj^)}. Then the volume of any finite polytope ?p bounded by these 
hyperplanes is a piecewise rational function of o,j, b, with integer coefficients. Indeed, the 
coordinates of the vertices of this polytope are (piecewise) rational functions of a,y, b; with 
integer coefficients according to Cramer's rule. Now triangulate the polytope *p into a set 
of X-dimensional simplices. The volume of any simplex with vertices (tiu, 021, . . . , t)Hi), 
1^1^ X+ 1, is equal to 

1 

Ol,3\f+l 
02,N+1 

Oxw+1 



1 

On 
D21 



1 

O12 
O22 
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Thus, the volume of *p is a polynomial in the coordinates of its vertices with coefficients 
which become integers after multiplication by N!. Various formulas and algorithms for com- 
puting the volumes of polytopes in Euclidean spaces of arbitrary dimensions are presented 
in the articles B6ll47ll48B9ll50l and references therein. 

A genuine proof of Conjecture [T]requires further studies beyond the present paper. 
8 Numerical experiments on the plane 

At previous stages of the project lllOlllll . we developed and implemented Fortran codes 
for computing all the hyperangular momenta J, K, A, and L (defined in Section [3]l and 
the energy terms T, Ta, Tp, T"', T\ T^, 7"="', 7'"', r"'=\ 7}, Tk, T^c, E""', EJ"", , E"\ 
iij", ™d £''=°"pi (defined in Sections |4] and |5]l for the physically interesting dimensions 
d = 2 and d = 3 of the ambient space. The input data for these codes are the entries of 
the d X n matrices Z and Z and the total mass M. For fixed dimension d, the number of 
operations required in the calculation grows linearly with the number of particles (N = n 
oiN = n + \),cf. RemarkU 

Using the codes prepared, we have verified Conjecture [T| (and, to some extent. Con- 
jecture |2j for systems of classical particles on the plane (d = 2) by numerical simulations. 
Within each of the two situations defined in Section |6] (particles with equal masses and par- 
ticles with random masses) and for each value of from 3 through 100 (i.e., for a total of 
98 values), £ = lO*" systems of A' particles on the Euclidean plane with the center-of-mass 
at the origin were chosen using a random number generator according to the procedures de- 
scribed in Section|6] A random point w uniformly distributed in the unit disc on centered 
at the origin was always generated as 

w = x^^^{cos(p, sin(p), 

where m and (p are independent random variables uniformly distributed in the intervals ^ 
X ^ I and 0^9^ 2k. For each system, we calculated all the energy terms T = I, Ta, Tp, 
r°\ P, T^, 7'="', r'"', Tj, Tk, T,^, E°''\ 4", and £™"p' defined in 

Sections |4]and|5] 

Similar simulations were performed in our previous paper [11] for J = 3 with 10^ sys- 
tems for each value of from 3 through 100 in each of the two situations. 

As one expects, for all the systems on M} in both the situations we found 7} = T™' 
(see Theoremlij and £'^"' = 0, i.e., fi™' = E""' [see ([35}]. For A' = 3 and all the systems in 
both the situations, we also found Tk = T"" (see Theorem |5]l and £" = 0, i.e., £)" = 
[see ([35]l]. 

The energies T^^^ , T^, and can be both positive and negative 1101 . We will denote 

their positive and negative "components" as 

Tl^' = max{r^' , 0) , TL^' = - min(r'='* , 0) , T'^' = Tl^' - Tl^' , 
T+ = max(rac , 0) , = - min(rac , 0) , T^c = T+ - T^ , 

4™P' = max(£™"Pi , 0) , = - min(£™"Pi , 0) , £™"p' = " • 

As was pointed out above, the quantities i?™' and Ef (and, consequently, E°"^, 
_^coupi^ and ii""'^') can be arbitrarily large for any fixed value T > 0. In our simulations, 
for any number A^ of particles, we encountered systems for which E°"^ > 5.42592 x 10"* 
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(for equal masses) or J5 5.1008 x 10"* (for random masses); recall that 2i™' = fj""' 
in our calculations. The maximal values of Z?""' (over all A') we observed turned out to be 
2.34534 X 10^ (for equal masses) and 2.0775 1 x 10^ (for random masses). Similarly, for any 
number A' of particles, there were systems for which E"^ > 5.41936 x lO'* (for equal masses) 
or £[" ^ 5. 10724 x 10^ (for random masses), and the maximal values of E'" over all A' were 
equal to 2.34534 x 10^ (for equal masses) and 2.07787 x 10^ (for random masses). Finally, 
for any number N of particles, systems occurred for which £''^™pi ^ —1.08453 x 10^ (for 
equal masses) or Zi'=°"pi ^ —1.0208 x 10^ (for random masses), and the minimal values of 
^coupi gygj. gjj ^ .^^gj-g _4,69068 X 10** (for equal masses) and -4.15538 x 10^ (for random 
masses). Moreover, in each of the two situations, the maximal values of Zi°"', Zi™, E^", 
and — ZJ"^""?' were attained at the same system — for which the "angle" between the spaces 
nieft(Z) and nright(Z) is very small (cf. C^). Since the quantities i'""', fij", and 
are unbounded, we did not examine their statistics. 

Remark 8 Of course, the large number of digits in the data above only reflects the particular 
set of numerical experiments. The same remark refers to similar data below where only 
several first digits are significant and informative. 

On the other hand, the energies TX, 7"', Tp, T', and do not exceed T = I while the 

terms 7'='" = 7>, T'"', r'=^ T^^', r'=^ 7>, 7;^ T+, T^, Ef, and "p' do not exceed 7™' 
according to Theorems |4] |5] and |8] For each of these quantities, we computed the mean 
values 

— 1 ^ 

^=£^^'' 

^1=1 

where Hi is the value of the energy H in question for the Ith system for the given N within 
the given situation (recall that £ = 10*), and the sample variances 



s^{H) = [H-Uf =H^-{Hf 

(see e.g. the manuals I51ll52ll53ll54l and references therein). Note that the sample variance is 
often defined as sl{H) = 2,s^{H)/{2,-\) I53II54I . The mathematical expectations ofs^{H) 
and sl(H) are equal to (£ - l)(Var//)/£ and VarZ/, respectively II51II52II53I I54I. where 
VarfZ = E [(/f - E//)^] = E(//-) - (E//)^ is the variance of H. However, for £ = 10*, the 
difference between the "biased" sample variance s^{H) and the "unbiased" sample variance 
sl(H) is of course negligible. 

The dependences of H on N for various terms H and both the mass distributions con- 
sidered are presented in Figs. [TJH] Along the abscissa axis on each of these figures, the 
"physical" distance x between the left end point corresponding to N = 3 and the point cor- 
responding to a given is proportional to j — . In such a coordinate frame, any depen- 
dence 

av + b ^, ^ 

y= — — , v=N-l, 
cV 

is represented by a straight line: 

I 1 2 _av + b _2a + b-2bx 

V ^ ~ 1 - 2jc ^ CV ~ 2c ■ 

In Figs. 12}^ we also show the fractions of systems for which T'^'^'^ < 0, Tac < 0, or £''^°"pi > 0. 
The "oscillations" on the corresponding dotted lines in Figs.|2]and|4]in the region of large 
N are due to some subtle shortcomings of the graphic system we used (gnuplot 4.0). 
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Fig. 1 Tlie mean values of tlie terms of tlie Smith decomposition 1141 and orthogonal decomposition l IlSl of 
the total kinetic energy T = 1 , as well as the mean values of the shape energy . The left panel corresponds 
to the equal mass case and the right one, to the random mass case. The physical space dimension d is equal 
to 2. 



For d = 2 and N ^ 3 (and, consequentiy, v = A' — 1 ^2 and CO = 2), the formulas 
(I27l)-ll3ni of Conjecture [T]take the form 

ETa = 1-^, ETp = :^, Er°^ = l-^, ET' = ^, (41) 

Er^ = 2^, (42) 

Er""' = — , Er"" = 1 - — , ET"''=0, (43) 
2v 2v 

m=^, ETk = ^, Er,c = ^, (44) 
E£2°"' = 0, E^" = 1 - ^ • (45) 

Of course, the equality E7''='" = E7> follows from the fact that 7'="' = Tj for d = 2, and the 
equality ES^"' = follows from the fact that generically £:|"' = for > <i + 1 . These 
two equalities hold for any distribution of the masses. Similarly, the equality EE™ = for 
V = 2 (i.e., for N = 3) follows from the fact that generically Z?™ = for N d +1 (again 
for any distribution of the masses). On the other hand, the formulas (I41b-( l43t indicate that 
ETp = ET^ = EJ'^"' for d = 2 and N ^ 3 in the situation of equal masses (although the 
hyperradial energy Tp , the shape energy , and the external energy T^^^ of an individual 
system on the plane are not related to each other at all). 

Thus, what we verified was the formulas (l41 M44t for the 10 terms T/^,Tp, 7"'°', T' , T^, 
jext^ r'"', r''\ Tk, T^c in the situation of equal masses f or 3 ^ A? sj 100, the formula (|45J 
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Fig. 2 The mean values of the terms of the projective partition ( 1161 of the rotational energy T"" for 7=1. 
The left panel corresponds to the equal mass case and the right one, to the random mass case. The physical 
space dimension d is equal to 2, so that T"^"' = 7j according to TheoreiTi|5] The mean values of T^_^'^ and 
yrcs indistinguishable on the scale of the figure, so that the inean values of the residual energy T"^^ itself 
are indistinguishable from zero. The dotted line labeled as "(7'^'='^ > 0) = {T'"^^ < 0)" presents the fraction of 
systems for which T"'" is positive or the fraction of systems for which T"'" is negative (both these fractions 
are indistinguishable from 1/2 on the scale of the figure). 



for the term E'" in the situation of equal masses for 4 ^ A' ^ 100, and the formula ^T'^'"^ = 
in the situation of random masses for 3 ^ 100 (see Conjecture |2]i — a total of 

11 X98 + 97 = 1175 

equalities of the form H = EH where the EH values are given by the formulas (I41b -( l45t . 
All these equalities were indeed obtained in our calculations. Namely, the minimal value, 
the maximal value, and the mean value of the difference \H — E//[ [over the 1 175 triples 

(the term H, N, equal/random masses) 

considered] were found to be 0, 4.81437 x 10""^, and 2.7892 x 10"^ respectively, and the 
minimal value, the maximal value, and the mean value of the ratio 

\h-eh\ . 

— ■ ' — £ = 10 

[.2(//)/£]'/^' 

turned out to be equal to 0, 2.41227, and 0.698384, respectively (so that indeed H = EH 
up to statistical errors Il51ll521l53ll54l . see Figs. [T}0}. Note also that the minimal value, 
the maximal value, and the mean value of the "weighted" difference 2v\H — EH\ were 
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Fig. 3 Tlie mean values of tlie terms of tiie hypersplierical partition 1171 of the rotational energy 7™' for 
r = 1. The left panel corresponds to the equal mass case and the right one. to the random mass case. The 
physical space dimension d is equal to 2, so that Tj = T"^"' according to Theorem [3] The dotted line labeled 
as "(Tac < 0)" presents the fraction of systems for which the angular coupling energy Tac is negative. 



0, 2.31975 X 10 ^, and 1.63148 x 10 ^, respectively, while the minimal value, the max- 
imal value, and the mean value of the sample variance s^{H) equaled 5.02904 x 10^^, 
8.45973 X 10"^ and 3.5151 x 10"^ respectively. 

We do not show the standard errors 8H = [s^{H) /£] = \Q-^s{H) of the mean values 
H 151,52 53 54 1 in Figs. since the numbers H + 5H and H — 5H are indistinguishable 
on the scale of the figures. 

The minimal value, the maximal value, and the mean value of the number of systems 
for which T"'* < (over the 98 values of N) were equal to 498853, 501308, and 499930 
(of £ = 10*), respectively, in the situation of equal masses and to 498691, 501 145, and 
499963, respectively, in the situation of random masses. Therefore, in both the situations, 
T'^"^ is negative for half of the systems up to statistical errors (see Fig.|2]i. 

We have thus confirmed Conjectures [T] and |2] for d = 2 and 3 ^ A' $J 100 (to be more 
precise, Conjecture|2]has been confirmed only in the situations of equal and random masses). 
Moreover, it turns out that s^{T^^) = 5^(7^^^) at any A' for both the mass distributions up to 
statistical errors. 

The equality EJ''^' = implies that E7^'=' = ETL^'' for both the mass distributions. Our 
numerical simulations suggest that for sufficiently large A' (say, for A' > 30) in the situation 
of equal masses, 

and this approximation is quite good. No similar formula exists for E,T^^ and ET'^^ in the 
situation of random masses. Recall that in our previous paper 1 11 1, we found that for spatial 
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Fig. 4 The mean values of the bounded terms of the singular value expansion (Ts) [and of its refinement J21H 
of the rotational energy 7'°' for T = 1. The left panel corresponds to the equal mass case and the right one, to 
the random mass case. The physical space dimension d is equal to 2. The dotted line labeled as "(£''^°"P' > 0)" 
presents the fraction of systems for which the tangent coupling energy £''^""P' is positive. 



systems of particles (d = 3), KT'^'^ = ET^^ w 5/(12v) for sufficiently large A' in the situation 
of equal masses. It would be interesting to figure out how the coefficients 5/16 and 5/12 
can be generalized to larger dimensions d. For instance, the conjectural formula 



8dv 



gives the desired values and for d = 2 and d = 3, respectively, and agrees with the 
fact that r"^** = for (i = 1 (independently of the masses of the particles, see Example [3ll. 

There are no analogs of the formulas (I4U - (I45) for the situation of random masses (ex- 
cept for the equalities ET'"'''^ = Eiij"' = 0, of course). However, for large A' (starting with 
N w 50, say), the following approximate equalities hold in our simulations in the situation 
of random masses: 

_ 1.88Af-4.04 _ 0.12Af + 2.04 

— 1.83A?-5.2 -J Q.llN + 'i.l — O.OSA'+l.lS 

2v ' 

_ _ 0.12iV + 2.05 -7- 1.7W-7.22 

yext — f T^mt ^ 

2v ' 2v ' 

0.03iV + 0.79 



yres _ yres 



2v 

_ 0.88A'-3.03 ^ 0.83A?-4.22 
7> « 



2v ' 2v 
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The coefficients a and b in these expressions H w {aN + b) / (Iv) are obtained for each mean 
energy H via minimizing the sum 

100 _ 

[2(A?-l)//(A?)-aA'-fc] . 

^=50 

In fact, we found b = 0.78 for Tp and & = 0.8 for T^. 

Since Tj = T'^"' for any A' (provided that d = 2) and Tk = T'"' for A'^ = 3, one concludes 
that = T'^'^^ for A' = 3 independently of the masses. Accordingly, in our simulations with 
A' = 3, Jac = and T^c is negative for half of the systems up to statistical errors for both the 
mass distributions. For larger values of A', one finds > 0, while the fraction of systems 
for which T^c is negative decreases rapidly as A^ grows for both the mass distributions (see 
Fig.O. However, for each N ^ 4, this fraction is greater in the situation of random masses. 
For particles with random masses, we encountered systems with a negative energy for 
any value of A'. On the other hand, for particles with equal masses for N = 32, N = 33, and 
37 ^ N ^ 100, the angular coupling energy J^c of ''H the 10^ systems we sampled turned 
out to be positive (for each value A' = 34, 35, and 36, we observed only one system with 
Tac < 0). In fact, systems with the center-of-mass at the origin and with a negative energy 

exist for any d ^ 2, N ^ 3 and any masses m\,m2,..., mf^, but the "relative measure" of 
the set of such systems is very small for large A' and mi = m2 = • • • = m^. To "construct" a 
system of A' > 3 particles with prescribed masses and with T^c < 0, one may choose a pair 
(Z, Z) corresponding to N = 3 and Tac < and then augment both the matrices Z and Z on 
the right by A^ — 3 zero columns (cf. 1 11]). 

For both the mass distributions and any value of A', the fraction of systems for which 
^jcoupi > is much less than 1 /2 and decreases fast as A' grows for both the mass distri- 
butions (see Fig.|4ll. For each A^, this fraction is greater in the situation of random masses. 
Nevertheless, even for A^ = 100 in the situation with equal masses, the tangent coupling 
energy £''=°"pi was positive for 32969 systems of 10''. 

In our numerical experiments for spatial systems of particles llll . the energies Tac and 
^coupi exhibited a similar behavior. 

In a complete analogy with the case d = 3 11 li . the 14 energy terms Tj<^,Tp, 7""', j', T^, 

7"="', J'"', Tk, Tac, T+, Tae , Ef, and in the present simulations for J = 2 can be 

divided into two classes, cf. Figs.[T}Q 

A) The seven terms Ta, T™', T'"', Tk, T^c, T.^, E™. For each of these energies H, 

the mean value H{N) increases as A' grows for both the mass distributions, and 

H(N) I , > H(N) I , , 3 ^ A? ^ 100. 

^ ' I equal masses ^ ' I random masses ' 

B) The seven terms Tp, T' , T^, r*^"', T^^, T^, Z?™"'''. For each of these energies H, 

the mean value H{N) decreases as A' grows for both the mass distributions, and 

— , — , (47) 

H(N) , <H(N) . , 3^A'^100. 

^ ' I equal masses ^ ' I random masses ' 

We do not consider here the terms r"^** and TI^^ because for both the mass distributions, 
yres _ Q jTCs _ jres ^j^y ^ statistical errors. We do not consider the term 7} 
either since Ti = T'^^^ for d = 2. 
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In some cases, the increase/decrease of H{N) as A' grows is slightly non-monotonous for 
large A' due to statistical errors. Moreover, in the situation of equal masses, as was pointed 
out above, turned out to be zero for all the systems sampled for N ^ 37 (in fact, was 
found to be very small for > 15). 

The exceptions to the rules (I46l>-ll47l) are as follows. First, E™ = for A' = 3 indepen- 
dently of the masses. Second, 



rr°t(3)| , =7™' (3) I , =i and accordingly 

^ ' I equal masses ^ ^ I random masses Z , 

— , — , , (48) 

^ ' I equal masses ^ ' I random masses 2 



up to Statistical errors. Third, 

T~(3)l =T~f3)l =0 

'^cv / I equal masses '^^^ (random masses 

up to statistical errors. Fourth, 

J7(N)\ , >7>(A')| . for 3^ A/ ^5 

S ^ ^1 equal masses S ^ ^1 random masses 

and 

I equal masses ^ I random masses 

within statistical errors. Fifth, 

JT(N)\ , >T^(N)\ , for 3 s; A^^ 5. 

-r \ y I equal masses i- \ / I random masses 



A and 

andom masses 



Sixth, 

J^(N)\ , >T^(N)\ 

'^^ ^ ' \ equal masses '^^ ^ ' I r 

ET^\m\ , >~^\n)\ , for 3^iV^4. 

+ ^ I equal masses + ^ ' I random masses 

In our previous paper Hill , being based on numerical simulations for spatial systems of 
particles (c? = 3) at 3 ^ A' ^ 4 and on the formulas (l33M34b concerning two-particle systems 
for any d, we conjectured that Er™' = (t/ - l)/<i and Er' = 1/t/ at all the values A' £/+ 1 
(and r = 1) for any mass distribution. The relations (I48t confirm this conjecture for d = 1, 
N = 3. 

The larger is the number A' of particles, the greater is the contribution of "kinematic 
rotations" to the total kinetic energy T of the system (see Remark [3]l and the smaller is the 
contribution of conventional rotations and changes in the singular values ^1,^2, ■ ■ ■ ,^a>- That 
is why the mean values of the energies J^, J"', T'"' (which include the contribution of 
"kinematic rotations"), and (which is connected with the kinematic angular momentum 
K) increase as A' grows for both the mass distributions whereas the mean values of the 
energies Tp, T' , , and 7"'''" = Tj decrease. It is however not clear why the mean values of 
the terms Tac, T^, and E2 increase as A^ grows while the mean values of the terms T^^, T.^, 
and ^'^""P' decrease. Why 

H(N)\ , >H(N)\ , for increasing terms H(N) (49) 

^ ^ I equal masses ^ ^ I random masses ^ ' ^ ^ 

and 

H(N)\ , <^(A')| ^ for decreasing terms H(N) (50) 

^ ^ I equal masses ^ ^ I random masses ^ ' ^ ^ 

is a complete mystery. 
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Similarly to the case J = 3 1111 . the sample variances s^{H) of each of the 15 energy 
terms Ta, Tp, 7™', T' , T^, 7"==", T'"', T"=\ r['=% Tk, T^c, T+, T^-, Ef, and "p' for d = 2 
decrease as the number of particles grows for both the mass distributions [recall that 
i2(7res-) ^ ^2(jre.s-, statistical errors]. For large A', this decrease is sometimes slightly 
non-monotonous due to statistical errors. Of course, s^{T.^) = in the situation of equal 
masses for > 37. As one expects, the sample variance of each of these terms in the situa- 
tion of random masses is larger than that in the situation of equal masses for the same value 
of N. For all the terms H, the sample variance s'^{H) gets very small in the situation of equal 
masses for large A'. 

The exceptions to these rules are as follows. First, in the situation of random masses, 
decreases starting with N = 4 (rather than with A' = 3): this quantity for A' = 4 is 
larger than for N = 3. Second, for both the mass distributions, s^{Tac) decreases starting with 
N = 4. Third, s^{T.^) decreases starting with A^ = 4 in the situation of equal masses and with 
A' = 5 in the situation of random masses. Fourth, ^^(Zij") decreases starting with A' = 4 in 
the situation of equal masses and with A' = 6 in the situation of random masses. Moreover, 
the variance of for A^ = 3 is of course zero independently of the masses. Apart from this, 
the inequality 

s^(H)\ , <s^iH)\ . (51) 

^ ^ I equal masses ^ ^ I random masses ^ ^ 

is violated in the following cases: 

for H = 7""^ and Tk at N = 3, 

for // = 4'' at A' = 4, 

for H = T^- and £•+ "p' at 3 ^ A' 4, 

for H = T^, r,e, and T+ at 3^Ni^5. 

As one sees, there is a strong correlation between the violation of the inequality i l51t and 
that of the inequalities flOll-dSOll. 

Our simulations confirm that the projective partition (I16t ensures a very effective sepa- 
ration between the conventional rotations and "kinematic rotations" l8l fTT]|15lll6l compared 
with the hyperspherical partition (I17K not to mention the singular value expansion dlSt . 
The mean absolute value | T'''^ \ of the residual energy decreases as N grows for both the 
mass distributions, whereas the mean absolute value \T^c\ of the angular coupling energy 
increases. For N = 3 one has T''"^ = T^^ independently of the masses, and consequently 
jres|('3^ = |rac|(3); these mean values are equal to 0.168796 in the situation of equal masses 
and to 0.146767 in the situation of random masses. For N ^ 4, one has \T""^\ < \T.^c\, and 
the larger A', the greater is this difference. At A' = 4 



|r'-'=H(4) . =0.136738 < I J'-'''^ 1(4) , =0.151223 < 

I I \ / I random masses i i v / | equal masses 

|7U(4)| , =0.214097 < ^(4) I , =0.250362, 

I ' ^ 'I random masses ' ' ^ -^1 equal masses ' 



while at A' = 100 



|r«H(100) , =0.00640506 < I 1(100) , =0.0381073 < 

I I ^ ^ I equal masses i i ^ •'I random masses 

170(100) I , =0.400477 <|70( 100) I , =0.489895. 

1 o.^ I \ / I random masses i i \ / \ equal masses 
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9 Conclusions 

The statistical studies of our previous paper lllll (devoted to systems in M.^} and those of 
the present paper (devoted to systems on M?) are formal in the sense that they do not take 
into account any interaction potentials between the particles. If one considers kinetic energy 
partitions for interacting particles with a certain potential energy U, then it is more natural 
to average various energy terms H at a fixed total energy T + U (averaging over a micro- 
canonical ensemble, see e.g. |55J) rather than at a fixed total kinetic energy T , cf. (5][6l[8l 
[T5III6I . Choosing potential energy hypersurfaces at random according to some distribution 
in an appropriate infinite dimensional functional space, one would probably obtain entirely 
different statistics of the energy components. In this case, it seems suitable to average the 
ratios H /T or ///(T + 11,) over the initial conditions, the potential U, and the time. It is also 
of interest to compute the mean values of the energy terms at fixed values of the total angu- 
lar momentum / (which is, by the way, customary in quantum mechanics, see e.g. 1561 ") or 
kinematic angular momentum K. 

There are many ways to generalize the energy partitions treated in the present work. 
One of them is pointed out in our previous paper 0111 and consists in defining the energy 
terms corresponding to the actions of arbitrary subgroups of the orthogonal groups 0{d) and 
O(A'). Another approach recently proposed by Marsden and coworkers 16J for J = 3 is called 
the hyperspherical mode analysis by the authors. The 2>N — 6 internal modes of an A'-atom 
system {N ^ 5) in are classified in f6l into three gyration-radius modes, three twisting 
modes, and 3A^ — 12 shearing modes. Most probably, Marsden's theory can be generalized 
to the case of an arbitrary dimension d. 

One of the main results of our previous paper fill and the present paper is that in the 
situation of equal masses, the mean values E// of various components H of the total kinetic 
energy T are expressed in terms of the dimension d of the physical space and the number 
N of particles in a very simple way. However, it is not clear at all whether the distributions 
of H are "simple" functions of d, N, and H, not to mention the joint distributions of several 
components. For instance, we have not attempted to find any expressions for the variances 
Var// of H [or, equivalently, for E(//^)] or, say, for the correlation coefficients II51II52||531 
l54l between the energy terms. We hope that such detailed statistical properties of the kinetic 
energy partitions of classical systems will be examined (both numerically and rigorously) in 
further research. 

The work of MBS was supported in part by a grant of the President of the Russia Fed- 
eration, project No. NSh-4850.2012.1. 
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